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Abstract: We provide an update of the overview of imprints of Tsallis nonextensive statistics seen in a multiparticle 

production processes. They reveal an ubiquitous presence of power law distributions of different variables 
characterized by the nonextensivity parameter q > 1. In nuclear collisions one additionally observes a in- 
dependence of the multiplicity fluctuations reflecting the finiteness of the hadronizing source. We present 
sum rules connecting parameters q obtained from an analysis of different observables, which allows us to 
combine different kinds of fluctuations seen in the data and analyze an ensemble in which the energy (E), 
temperature (T) and multiplicity (N) can all fluctuate. This results in a generalization of the so called 
Lindhard's thermodynamic uncertainty relation. Finally, based on the example of nucleus-nucleus collisions 
(treated as a quasi-superposition of nucleon-nucleon collisions) we demonstrate that, for the standard 
Tsallis entropy with degree of nonextensivity q < 1, the corresponding standard Tsallis distribution is 
described by q' = 2 — q > 1. 
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1. Introduction 

Multiparticle production processes are the main source of our knowledge of the properties of matter formed in 
extreme conditions. In fact, they are sometimes regarded as an illustration, in laboratory conditions, of the 
state apparently existing just after the Big Bang. From the very beginning for their investigation it was obvious 
that the best way of their description and understanding is the approach of statistical mechanics 1 . For a decade 
now we advocate (see review [2] for references) that the quality of data is high enough to observe the deviations 
from the usual statistical approach based on a Boltzman-Gibbs ensemble towards a more refined one based on 
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1 See, for example [1] and references therein. 



The imprints of supcrstatistics in multiparticle production processes 



its generalization, for example given by Tsallis statistics [3-5] 2 . From the very beginning, we have identified 
the q parameter occurring in Tsallis statistics with the measure of some intrinsic fluctuations present in the 
system [6, 7]. They are natural in systems where the heat bath is not homogeneous. In short: if the scale 



parameter T ("temperature") in the usual Boltzman distribution, f{E) oc exp y—E/Tj, fluctuates according to 
a gamma distribution, g(l/f) oc |t/ — l)rj j exp |— T/ ^{q — 1)t| j, with fluctuations given by 

g = l + Var(l/f)/{l/f) 2 , as result one gets the power-like Tsallis distribution: 



2 r'^J 



1-(1-5)| 



(1) 



This was quickly recognized as an example of so called superstatistics [8, 9] and, in what follows, we shall 
understand it in this way, i.e., we shall mainly concentrate on the notion of fluctuations represented by the 
parameter q. 

The conclusion of [6, 7] was reinforced by a more refined analysis in [10, 11] and generalized in [2] by additionally 
allowing for the heat bath to exchange energy with its surroundings. One then gets the same Tsallis distribution 
as before, but with a g-dependent effective temperature, T e ff = T+ (q — 1)T», replacing T in Eq. (1). Here T* is a 
new parameter depending on the transport properties of the space surrounding the emission region. In [2] , where 
this concept was first introduced for heavy ion collisions, T, = 4>/{Dcpp), cf> is the energy transfer taking place 
between the source and its surroundings 3 . This quantity is supposed to model the possible transfer of energy from 
the central region of nucleus-nucleus interaction towards the spectators not participating in the collision (when 

r. < o) 4 . 

For nuclear collisions the possible transfer of energy between the region of interaction and the spectators depends 
on its size (or on centrality of the reaction). Accounting for it thus explains automatically the observed q- 
dependence of the multiplicity fluctuations measured in nuclear collisions on the centrality of collision [12] (cf. 
also [13]). Assuming that the size of the thermal system produced in heavy ion collisions is proportional to the 
number of nucleons participating in collision, Np, one gets that (CV is heat capacity under constant volume) 

which nicely fits the data [12]. 

Since the review [2], Tsallis distributions have been used more widely, for example, they were successfully applied 
to the analysis of the recent PHENIX [15] and LHC CMS [16, 17] data (cf. also compilation [13]). In what 



2 For an updated bibliography on this subject see http://tsallis.cat.cbpf.br/biblio.htm. 

3 It can be connected with viscosity r\ by <f> = rjf(u), where /(it) = (dui/dxk + du^/dxi) 2 (u being velocity); D, 
cp and p are, respectively, the strength of the temperature fluctuations, the specific heat under constant pressure 
and the density. 

4 It is interesting to note that in cosmic ray physics it can be argued [14] that the corresponding T* > and 
describes the transfer of energy from the surroundings to cosmic ray particles accelerated in outer space. 
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follows we present our recent results in this field in more detail: a derivation of g-sum rules unifying different 
fluctuations - in Section 2, a generalization of thermodynamic uncertainty relations and their applications - in 
Section 3. Section 4 contains our present result, namely, using as an example the experimental data on nucleus- 
nucleus collisions (treated as quasi-superposition of nucleon-nucleon collisions), we demonstrate the interrelation 
between nonextensivity parameters obtained from Tsallis entropy and Tsallis distributions. Section 5 summarizes 
our work. 

2. g-sum rules 

As argued above, the parameter q provides a useful measure of intrinsic (nonstatistical) fluctuations in the system 
[6]. In multiparticle production processes the particles produced are characterized by their positions in phase 
space, in particular the measured momentum p is decomposed into longitudinal component (along the direction 
of colliding particles), p^ — mysinhj/, and transverse component, px, with p — \J\p 2 \ = \Jp\ + p\ (here 
m,T = \/ vn? + p\ denotes the so called " transverse mass" of the particle and y = | In %^ P p L L its " rapidity" ; E is 
the energy of the particle and m is its mass) . Because in most cases data are presented in the form of distributions 
in rapidity y (i.e., they are integrated overpr), dN/dy, and as distributions inpr (i.e., they are integrated over y), 
dN/dpx, one is confronted with two different fluctuations: in longitudinal phase space, characterized by q = q^, 
and in transverse phase space characterized by q = qr- ft turns out that the strengths of both fluctuations 
measured by q is different, whereas §l — 1 ~ 0.1 — 0.3 and grows with energy of collision (measured mainly in 
pp and pp collisions), transverse fluctuations are much weaker, qr — 1 ~ 0.01 — 0.1, vary slowly with energy and 
depend slightly on whether one observes elementary collisions or collisions between nuclei [19-22]. 
There is another source of knowledge concerning fluctuations, namely observed multiplicity distributions, P(N). It 
turns out that temperature fluctuations in the form of a gamma distribution leading to Eq. (1) result in substantial 
broadening of the corresponding multiplicity distributions. This changes from poissonian form characteristic for 
exponential distributions, P(N) = N exp (— iV) /N\ (where N = E/T), to negative binomial form (NB) for 
Tsallis distributions, Eq. (1) (cf., [23], for details) 5 , 



r(iv + i)r(fc) ( 1 + i^>y 



The nonextensivity parameter q reflects here fluctuations in the whole of phase space and can be (and usually is), 
different from the previously obtained qL and qr- In [19, 20] it was proposed that because q — 1 = a 2 (T)/T 2 (i.e., 
it is given by fluctuations of total temperature T), then assuming that a 2 (T) = o 2 (Tl) + (J 2 (Tt), the resulting 



5 Notice that in the limiting cases of q — > 1 one has k — > oo and (3) becomes a poissonian distribution, whereas 
for q — > 2 on has k — > 1 and (3) becomes a geometrical distribution. 
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values of q should not be too different from 

q = + _ Tj + n + L (4) 

Therefore for the observed dominance of longitudinal (partition) temperature over the transverse one, Tl 3> Tt, 
one should expect that q ~ which is indeed observed [19, 20]. This is the first sum rule observed for parameters 
q obtained from different measurements. 

In cases where other variables in addition to T also fluctuate (usually fluctuations of temperature are deduced 
either from data averaged over other fluctuations or from data accounting also for fluctuations of other variables) 
one should refine the experimentally evaluated q. In this case, when extracting q from distributions of dN/dy, 
one finds that (cf., [24] for details) 

def Var(T) _ Var(z) Var(m T ) _ m T , . 

q ~ <T> 2 ~ (m T y ; Z ~^T- ( ) 

This is the second sum rule for parameters q obtained from different measurements. It connects the total q, 
which can be obtained from an analysis of the NB form of the measured multiplicity distributions, P(N), with 
qL — 1 = Var(z)/ (z) 2 , obtained from fitting rapidity distributions and Var (wit) /{wit) 2 obtained from data on 
transverse mass distributions. When extracting q from distributions of dN/drriT, we proceed analogously with 
2 = coshy/T. 

3. Generalized thermodynamic uncertainty relations 

We shall continue the above discussion introducing the notion of thermodynamic uncertainty relations and propos- 
ing their generalization with the help of nonextensive statistics [25]. They were discussed in [26] where it was 
suggested that the temperature T and energy U could be regarded as being complementary in the same way as 
energy and time are in quantum mechanics. A simple dimensional analysis suggests that AU A/3 > k, where 
P — 1/T and k is Boltzmann's constant. Isolation (U definite) and contact with a heat bath (T definite) are then 
the two extreme cases of such complementarity. This is known as Lindhard's uncertainty relation between the 
fluctuations of U and T [27]: 

LOu + u)' T = where w 2 = Var{x)/(x) 2 . (6) 

This idea is still disputable [28-30] , nevertheless we can treat these increments as a measure of fluctuations of the 
corresponding physical quantities. This allows us to analyze an ensemble in which the energy (U), temperature 
(T) and multiplicity (N), can all fluctuate and thus to express these fluctuations by the corresponding parameters 
q. In this way, using generalized thermodynamics (based on nonextensive statistics) one gets the following relation 
[25] 

| W iV — t^t| = 0Ju + oj? — 2pLouu) T = (wc; — u) T ) 2 + 2ujjju>t{1 — p) = \q — 1|, (7) 
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Figure 1. (Color online) Energy dependencies of the parameters q obtained from different observablcs. Triangles: q obtained from 
an analysis of rapidity distributions [18, 31]. Squares: q obtained from multiplicity distributions P(N) [32, 33] (fitted 
by q — 1 4- 1/fc with 1/k — —0.104 + 0.0291ns). Circles: q obtained from different analysis of transverse momenta 
distribution /(pt). Data points in this case come, respectively, from the compilation of p + P data (full symbols) [34], 
from CMS data (half filled circles at high energies) [16, 17] and from NA49 data from Pb+Pb collisions (half filled 
circles at low energies) [35—37]. The full line comes from our recalculation using Eq.(7). 

where p — p(U,T) g [— 1, 1] is the correlation coefficient between U and T. This generalizes Linhard's thermody- 
namic uncertainty relation. 

The observed systematics in energy dependence of the parameter q deduced from presently available data is shown 
in Fig. 1. From measurements of different observables one observes that for high enough energies q > 1 (for low 
energies conservation laws are important and one can encounter q < 1 situations) and that values of q found from 
different observables are different. The later is caused either by technical (methodical) problems or else from some 
physical cause. The former arises when, for example, fluctuations of temperature are deduced either from data 
averaged over other fluctuations or from more refined data also accounting for fluctuations of other variables (as 
in [24], see Eq. (5)). The latter case is connected with the fact that the observed q's were obtained in different 
parts of phase space (or both). In this case one gets an uncertainty relation (7) with the help of which one can 
connect fluctuations observed in different parts of phase space. For example, one can recalculate q obtained from 
P(N) (dashed line in Fig. I) to q which can be evaluated from f(pr) (full line in Fig. I), see [25] for details 6 . 

6 A comment is necessary when looking at results at Fig. 1 obtained from f(y). Namely, as observed in [13, 24], 
it turns out that in the fitting procedure parameters T and q are strongly correlated. This is why q values evaluated 
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4. Nonadditivity in nuclear collisions (on the q duality in nonex- 
tensive statistics) 

So far we were using Tsallis statistics without really resorting to Tsallis entropy, i.e., we treated it as a kind of 
superstatistics [8, 9]. However, closer inspection of both approaches reveals that corresponding nonextensivity 
parameters (say q and q' , respectively) are not the same, in fact one encounters a sort of duality, like q — 2 — q 
discussed, for example, in [10, 11, 38]. We shall now address this problem in more detail on the example of 
nonadditivity observed in nuclear collisions where, as we shall see, both types of q can be discussed and (in 
principle) compared at the same time 7 . 

One of the phenomenological approaches used to describe these collisions is based on superposition models in 
which the main ingredients are nucleons which have interacted at least once [42]. In this case, when sources are 
identical and independent of each other, the total (N) and the mean ({N)) multiplicities are supposed to be given 
by, 

N = ^2m, and (N) = (u)(m), (8) 
!=i 

where v denotes the number of sources and n; the multiplicity of secondaries from the i th source. Albeit at present 
nuclear collisions are mostly described by different kinds of statistical models [1] , which automatically account for 
possible collective effects, nevertheless a surprisingly large amount of data can still be described by assuming the 
above superposition of independent nucleon-nucleon collisions (possibly slightly modified) as the main mechanism 
of production of secondaries and the question of the range of its validity is a legitimate one [43, 44]. 
Using the notion of entropy, and considering v independent systems for which the corresponding individual 
probabilities are combined as 

i> 

P M( Ii ,..., I( „)) = []p1 1) W, (9) 

fc=i 

and assuming that all p q (xk) are the same for all k (i.e., their corresponding entropies Sq 1 ^ are equal), one finds 
that 8 



k 



-J (v-k)\kl 



in different analysis of rapidity distributions [18, 31] differ slightly from presented here (roughly speaking, they 
give q values comparable or something higher that one obtained from multiplicity distribution). 

7 Similar duality occurs in nonextensive treatment of fermions for which the particle-hole correspondence, 
n q (E,T, /i) = 1 — n2- q (—E,T,—jj,) (where /i is the chemical potential), must be preserved by the q-Fermi dis- 
tributions [39, 40]- However, here we deal with different problem, namely that parameter q in entropy S q differs 
from parameter q' in probability distribution f q i and that q — 2 — q . 

8 Notice that In [l + (1 - g)4' y) ] = i/ln [l + (1 - q)S q 1] ^ and S { q ] q -^l v ■ Sf> . For q < 1 one has S { q ] /v oo, 

i.e., entropy Sq is nonextensive. For q > 1 one has S q ^ > only for q < 1 + l/S, 1 ' and S q jv u zH^ o, i.e., 
entropy is extensive, < Sq jv < S q r K 
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In the following we put v = Nw/2 = Np (Nw is the number of wounded nucleons and Np is the number of 
participants from a projectile). Assuming naively that the total entropy is proportional to the mean number of 
produced particles, 

S = a{N), (11) 
one obtains the following relation between mean multiplicities in AA and NN collisions, 

a(N) AA = [1 + (1-9)^^-1. (12) 

At this point we stress the following important observation, so far not discussed in detail. Namely, because (as 
shown in [12]), {N)aa increases nonlinearly with Np and {N)aa > Np ■ {N) pp , the nonextensivity parameter 
obtained here from considering the corresponding entropies must be smaller than unity, q < 1. On the other 
hand, all estimations od the nonextensivity parameter (let us denote it by q') discussed before lead to q' > 1. 
This is an apparent q duality in nonextensive statistics, on which we shall concentrate in more detail. 
Start with the obvious remark that, strictly speaking, relation (12) is not exactly correct for S q . In what follows 
we denote entropy on the level of particle production by s (and the corresponding nonextensivity parameter by 
q), whereas the corresponding entropies and nonextensivity parameter on the level of NN collisions by S and q. 
From Eq. (10) we have that for N particles 



[l + (!-«) sf] 



" '' I ^ 9 >,1 „ (1) , r ,-. Q N 

s\ = -t ^ > N ■ sy = aN, (13) 



where s^ = a is the entropy of a single particle. In a A + A collision with v nucleons participating, Eq. (10) 
results in 

[l + (l 



= ± ^— ^ , (14) 



where Sq is the entropy of a single nucleon. 

Denoting multiplicity in single N + N collision by n, the respective entropy is Sq = S^' 



(AT) 
9 



1 + (1 — q) s|^j — l|/(l — <j), whereas entropy in A + A collision for N produced particles is S 
I [l + (1 - q) s^'J N - l| /(l - q). This means that 

Sf>=S<?\ (15) 

Notice that parameters q and q are usually not identical. Moreover, from relation (2) one gets that for iViV 
collisions (where Np — A) q — 1. On the other hand, for q = q Eq. (15) corresponds to the situation encountered 
in superpositions, as in this case one gets that 

l + (l-g)4 1) ] JV = + Q)4 1> Y L/ or N = nv. (16) 
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Fi gure 2. (Color online) Energy dependence of the charged multiplicity for nucleus-nucleus collisions divided by the superposition of 
multiplicities from proton-proton collisions (cf. Eq. (19)). Experimental data on multiplicity arc taken from compilation 
[41]. 
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Fi glire 3. (Color online) Energy (-\/s) dependence of the parameter ci(s) obtained from the analysis of available data 
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Consider now the general case and denote 



Ci = 1 + (1 - q) 4 



(i). 



C2 



l-q 
1-q 



(17) 



These quantities are not independent because: 



C2Cx + 1 — C2 = (C2C a + 1 — C 2 ) . 



(18) 



From relation (18) one gets that 



N 



v ■ n vn ■ In a 



In 



(c 2 c? + 1 - c 2 y - (1 - oa) 



f-'L' 



(19) 



which for iV = (Naa), n = {N pp ) and ^ = iVp is presented in Fig. 2 for different reactions. As seen there one 
can describe experimental data by using c 2 = 1.7 and with ci depending on energy yfs according to ci(s) = 



1.0006 - 0.036s" 



as seen in Fig. 3. Notice that for energies y's > 7 GeV one has ci > 1. This means that 



q < 1 and (because c 2 > 0) also q < 1. 

Now look at this problem from the view point of Tsallis entropy, 



1 



dxf q (x) - 1 



(20) 



To get from it the probability density function fix), one either optimizes it with constrains 



j dxf(x) = 1; J dxxf q (x) = 



and obtains [45] 



f(x) = (2-q)[l-(l-q)x]*=*i < x < oo; 1 < q < 3/2, 
or else one uses as constrains 



J dxf(x) — 1; ^ dxxf(x) = (a;) 



and obtains [45] 



/(*) = 



< a; < oo; 1/2 < g < 1. 



(21) 



(22) 



(23) 



(24) 



[1 + (1 - q)x] i-a 

Notice now that only (22) is the same as distribution obtained in superstatistics and used above, cf., Eq. (1). The 
second distribution, Eq. (24), which seems to be more natural from the point of view of a physical interpretation 
of the constraint used, becomes the first one if expressed in q' given by 
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namely, in this case one has 

f(x) = (2 - q') [1 - (1 - q')x] ^ . (26) 

We show here, cf. Fig. 1, that using a Tsallis distribution in the form of Eq. (26), one gets q' > 1. On the other 
hand, non additivity in the superposition model described using the notion of entropy clearly requires q < 1, cf. 
Figs. 2 and 3. This means that q is not the same as q. The conclusion one can derive from these considerations 
is that the second way of deriving f(x), which uses a linear condition, cf. Eq. (24), is the correct one and that 
q' in distribution is not the same as q in entropy. The problem is that, whereas from distributions one can easily 
deduce a numerical value of q' , this is not the case when one uses entropy. There are too many variables to play 
with (cf., considerations using the superposition model as above). For example, in the definition of ci in Eq. (17), 
one has the , which is not known a priori. The only thing one can get in this case is that q < 1. We cannot 
therefore check numerically that relation (25) really holds. But, if one agrees that the Tsallis distribution comes 
from Tsallis entropy, we have only two options: either q' = q or q — 1 = 1 — q. Our conclusion presented here, 
that q' > 1 and q < 1, therefore supports the second option, i.e., Eq. (25). 

5. Summary 

To summarize, Tsallis statistics is fruitful because in a very economical way (with only one new parameter 
q) it describes the power-like behavior of different observables. This parameter, for q > 1 considered here, is 
given fully by nonstatistical fluctuations present in the system and visible as fluctuations of the scale parameter 
in superstatistics. It also allows (via specific sum rules or through a generalized thermodynamic uncertainty 
relation) to connect fluctuations of different observables or observed in different parts of phase space. Finally, 
when considering a superposition scenario, for example, in the scattering of nuclei, the relation q' — 1 = 1 — q 
seems to be observed (with q' occurring in the Tsallis distribution and q in Tsallis entropy). 

This final observation needs some more attention. The probability density function (PDF) is commonly evaluated 
by the Maximum Entropy Method (MEM) for Tsallis entropy with some constraints [46] 9 . At the moment, there 
are four possible MEMs discussed at length in [47] using two kinds of definition for an expectation value of 
physical quantities: the normal average (23) and the q-average (21) (with normal, as here, or the so-called escort 
PDFs [48-50]). Various arguments have been given justifying the g-average [51-53]. Recently, however, it has 
been pointed out that, for a small change of the PDF, thermodynamic averages obtained by the g-averages are 
unstable, whereas those obtained by the normal average are stable [54, 55]. On the other hand, it is claimed [56] 

9 Notice that Tsallis entropy is a monotonic function of the Renyi entropy, S q = ln 9 [exp (R q )], and that both lead 
to the same equilibrium statistics of particles (with coinciding maxima in equilibrium for similar constraints on 
the expectation value). 



in 
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that for the escort PDF, the Tsallis entropy and thermodynamical averages are robust. This means that this issue 
on the stability (robustness) of thermodynamical averages as well as the Tsallis entropy is still controversial [57] . 
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